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COOLING OF YOUNG STARS GROWING BY DISK ACCRETION. 
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ABSTRACT 

In the initial formation stages young stars must acquire a significant fraction of their mass by accretion 
from a circumstcllar disk that forms in the center of a collapsing protostellar cloud. Throughout this 
period mass accretion rates through the disk can reach 1CP 6 — 1CP 5 M Q yr _1 leading to substantial 
energy release in the vicinity of stellar surface. We study the impact of irradiation of the stellar surface 
produced by the hot inner disk on properties of accreting fully convective low-mass stars, and also look 
at objects such as young brown dwarfs and giant planets. At high accretion rates irradiation raises 
the surface temperature of the equatorial region above the photospheric temperature To that a star 
would have in the absence of accretion. The high-latitude (polar) parts of the stellar surface, where 
disk irradiation is weak, preserve their temperature at the level of To. In strongly irradiated regions 
an almost isothermal outer radiative zone forms on top of the fully convective interior, leading to 
the suppression of the local internal cooling flux derived from stellar contraction (similar suppression 
occurs in irradiated "hot Jupiters"). Properties of this radiative zone likely determine the amount of 
thermal energy that gets advected into the convective interior of the star. Total intrinsic luminosity 
integrated over the whole stellar surface is reduced compared to the non-accreting case, by up to a 
factor of several in some systems (young brown dwarfs, stars in quasar disks, forming giants planets), 
potentially leading to the retardation of stellar contraction. Stars and brown dwarfs irradiated by their 
disks tend to lose energy predominantly through their cool polar regions while young giant planets 
accreting through the disk cool through their whole surface. 

Subject headings: planets and satellites: formation — solar system: formation 
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1. INTRODUCTION. 

Our understanding of advanced stages of star forma- 
tion (T Tauri and later phases) has been significantly im- 
proved with the advent of infrared, submillimctcr, and 
high-resolution optical observatories such as HST and 
Spitzer. At the same time a great deal of uncertainty 
still remains regarding the earliest, so-called Class and 
Class I, stages of the star formation process. In the con- 
ventional nomenclature, Class stars are the protostellar 
cloud cores in the very beginning of their collapse, while 
Class I are the protostars embedded within an envelope 
of circumstellar material that is infalling, accumulating 
in the centrifugally supported disk, and being accreted 
by the protostars. 

At present, our knowledge exhibits a significant gap 
when it comes to describing the actual buildup of the 
stellar mass, from = Mq in the Class phase to 
M* ~ 1 Mq in the end of Class I phase. From the obser- 
vational point of view the major reasons for this are (1) 
the heavy obscuration provided by the increased densi- 
ties in the central part of the infalling protostellar core 
and the molecular cloud as a whole and (2) the difficulty 
in deriving the spectra of the central objects, namely 
distinguishing between the intrinsic protostellar and ac- 
cretion luminosities. At the same time, our ignorance 
concerns not only the history of mass accumulation by 
the protostars. Thermodynamical state of the accumu- 
lated gas is also an important ingredient of the picture. 
Stars that form out of material with high entropy, in par- 
ticular that processed through the accretion shock, tend 
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to have large sizes, while objects formed out of the lower 
entropy gas should be more compact. At the moment the 
uncertainty in the initial thermodynamical state of pro- 
tostellar objects precludes us from getting a good han- 
dle on the evolutionary tracks of the fully assembled (in 
terms of mass) protostars in the first 1 — 10 Myrs after 
their formation. Beyond about 10 Myrs, when the ini- 
tial conditions become largely forgotten, the evolution 
tracks calculated under different assumptions about the 
initial conditions typically converge (Baraffe et al. 2002). 
However, prior to this stage there are significant discrep- 
ancies between the results of different groups, and the 
uncertainty in the initial conditions for such calculations 
is the prime suspect for the difference. 

It is generally accepted that the conservation of 
the angular momentum in the collapsing protostellar 
cloud results in accumulation of the collapsed gas in a 
rotationally-supported disk in the cloud center. Only a 
small fraction of the cloud mass has low enough angular 
momentum to collapse directly into the protostellar core. 
The majority of stellar mass is most likely accumulated 
by accretion from the disk. According to observations, 
Class I stars acquire most of their mass on timescale of 
several 10 5 yrs which implies that disks around these 
objects must exhibit time averaged mass accretion rates 
of M ~ 10~ 6 — 10 -5 M Q yr _1 . Accretion luminosity re- 
leased in processing such large mass flux through the disk 
exceeds the intrinsic luminosity of the protostar. This 
immediately raises an issue of the possible non-trivial 
radiative coupling between the protostar and its circum- 
stellar disk. 

Effects of disk accretion on structure of young stars 
have been investigated by Mercer-Smith et al. (1984), 
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Fig. 1. — Schematic representation of stellar illumination by the 
disk (slant hashed). Filled region at the point where disk joins the 
star marks the boundary layer where intense energy dissipation 
takes place. Part of radiation emitted by the disk (arrows) gets 
intercepted by the star which heating it to temperature Ti rr , higher 
than the temperature To that a star would have in the absence of 
irradiation. Photospheric temperature is preserved at the level 
of To only in the polar regions of the star (marked with dashed 
lines) where disk illumination is weak. An external radiative zone 
(horizontally hashed) forms in the strongly irradiated parts of the 
stellar surface. 



Palla & Stahler (1992), Siess & Forestini (1996), Hart- 
mann et al. (1997), Siess et al. (1997, 1999). Some of 
these authors studies how the heat advected into the 
star with the freshly accreted material affects protostellar 
properties. However, none of these investigations looked 
at the effect of heat deposited at the stellar surface by 
radiation originating in the inner parts of the circumstel- 
lar disk, where most of the accretion energy is released 
(see Figure Q] for a schematic representation). Given 
that accretion luminosity may easily exceed the intrinsic 
stellar luminosity (luminosity derived from gravitational 
contraction, cooling and, possibly, deuterium burning in 
stellar interior), omission of this effect may not be justi- 
fied in many cases. 

In this paper we investigate stellar irradiation by the 
circumstellar disk and address the importance of this 
effect in determining the intrinsic luminosity of young 
stars. We calculate the spatial distribution of the disk 
flux on the stellar surface and determine when irradi- 
ation is important in ^2] The effect of irradiation on 
stellar cooling is investigated locally in $3] and globally 
in 21 Finally, in Sj5] we discuss the applications of this 
study to some real objects and its possible limitations. 

2. TEMPERATURE DISTRIBUTION DUE TO DISK IRRADIATION. 

We start by calculating the distribution on the stellar 
surface of the radiative flux Fi rr produced by the disk. 
We consider an axisymmetric geometrically thin disk ac- 
creting onto a star with radius i?* and mass M*. Flux 
F^r intercepted by the star is a function of 9 - the angle 
between the normal to the stellar surface and the normal 
to the disk (coincident with the polar axis of the star, 
assuming that disk lies in the stellar equatorial plane). 
Polar regions of the star are exposed to the radiation of 
only the distant, cool parts of the disk, while the equa- 
torial regions have a direct view to the innermost parts 



of the disk where most of the energy is dissipated. One 
can easily show that a disk extending all the way to the 
stellar surface gives rise to irradiation flux Fi rr {9) given 
by (Adams & Shu 1986; Popham 1997) 



Fi rr {0) 2 
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F d (R)RdR 
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{R 2 +B$- 2R+R sin 9 cos < 
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where R is the cylindrical radius, cos0 c = ii*/(r sin 6), 
Rin — R*/ cos#, and Fd(R) is the energy radiated by 
the unit surface area of the disk per unit of time. In 
Appendix A we demonstrate that this expression can be 
reduced to a one-dimensional integral which is easier to 
analyze than equation (fTj). 

To find the explicit dependence of Fi rr on 9 one needs 
to know Fd(R) which is determined by the viscous dissi- 
pation in the disk. Studies of steady-state thin accretion 
disks have generally found that 
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where M is a mass accretion rate and the function f(R), 
embodying the details of the disk emissivity near the 
stellar surface, behaves as / — > 1 when R ^S> R+- With 
Fd given by @ one finds 



, , GM+M , s 

F irr {e) = - 1 jt-g{0), 



(3) 



where the dimcnsionless function g{9) is given by equa- 
tion (fA3j) . 

A standard disk with zero torque at the stellar sur- 
face (situation appropriate for accretion onto black holes) 
has (Shakura & Sunyaev 1973) f(R) = 1 - (R±/R) 1/2 . 
The total viscous dissipation in such a disk is Ed = 
(l/2)GM ir M/R+ and the gas at the inner edge of the 
disk rotates at the local Keplerian velocity. This is in- 
appropriate in our case since the gas speed has to match 
the velocity of the stellar surface at R = R* (for sim- 
plicity assumed to be zero in our case). As a result 
a boundary layer must form near the stellar surface in 
which the azimuthal velocity of the gas is lowered by 
the viscous torque from the local Keplerian value to the 
stellar rotation speed. Viscous dissipation dramatically 
increases gas temperature in this layer creating an addi- 
tional source of radiative flux very close to the stellar sur- 
face. Irradiation by the boundary layer emission boosts 
up the stellar surface temperature in a narrow belt at the 
equator (with the thickness in ^-direction comparable to 
the thickness of the boundary layer) above that expected 
from the irradiation by the more distant parts of the disk, 
outside of the boundary layer. Thus, the existence of the 
boundary layer significantly modifies disk structure and 
emissivity near the stellar surface (Popham et al. 1993; 
Popham & Narayan 1995) complicating the calculation 
of/(i?). 

Fortunately, it will be shown later in 21 that cooling of 
irradiated stars depends only weakly on the behavior of 
f(R) at R ~ R+ and is thus relatively insensitive to the 
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Fig. 2. — Irradiation flux absorbed by the stellar surface in 
units of GM+M /R3 as a function of 6. Solid curve corresponds 
to Fd(R) oc il — 3 , dashed curve corresponds to Fd(R) oc R~ 3 [1 — 
(jR*/-R) 1/2 ], while the dotted curve represents asymptotic behavior 



structure of the boundary layer. For the mass accretion 
rates considered in this work (10~ 6 — 10 -5 M Q yr _1 ) the 
geometric thickness of the boundary layer is rather small, 
< 0.15 R* (Popham et al. 1993), so that the fraction of 
the stellar surface covered by the boundary layer and af- 
fected by the energy dissipation in it is rather small. For 
this reason we will further assume for simplicity that 3 
/(i?) w 1. In this case Fd keeps increasing all the way to 
the stellar surface [unlike the zero inner torque case in 
which Fd(R±) — * 0] thus roughly mimicking the contribu- 
tion of the boundary layer to the disk flux. We plot the 
behavior of function g{9) in Figure [2] for both f(R) = 1 
and f(R) = 1 - (R*/R) 1 ' 2 . 

Irrespective of the complications related to the exis- 
tence of the boundary layer one can derive useful results 
for Fi rr {9) in two asymptotic regimes. In particular, us- 
ing equation (jAlj) one finds that as 9 — > tt/2 



F 

J- 7 7 



F d (R*) 



fJ 



3 

16tt' 



(4) 



a result that is easy to understand since any point at 
the stellar equator receives disk radiation with uniform 
temperature corresponding to local disk flux Fd(R*) from 
7r steradian and reemits it into 2ir steradian. In the other 
limit of 6 — > one finds from equation (|A1[) that 



g » h sin 5 9, 



(5) 



where constant 1\ is given by equation (|A4[) . One can see 
from Figure [2] that approximation (J5]) works quite well 

3 Such assumption results in E d = (3/2)GM it M/R ir , larger than 
what can be provided by the change of the potential energy of disk 
material, but this inconsistency is not going to strongly affect our 
results. 



(better than 22% accuracy) for 9 < 0.5. This asymptotic 
behavior is insensitive to the details of the disk emissivity 
at R ~ R± since polar regions are irradiated only by parts 
of the disk at R > R± where f(R) w 1. 

Let us denote To and Lq = \-kR\uTq the temperature 
and luminosity which a star with mass M+ and radius i?* 
would have in the absence of irradiation (a is a Stephan- 
Boltzmann constant). To characterize the importance of 
irradiation we introduce irradiation parameter A: 



A = 



GM+M 



1.6 MtM-gR^T, 



3.5 > 



(6) 



where T n = T /10 n K, R n = J^/10™ cm, M x = M+/M G , 
and M n = M/(10™ M yr" 1 ). By construction, A is 
roughly the ratio of accretion luminosity and the stellar 
luminosity Lq in the absence of irradiation. 

When irradiation is allowed for the photospheric tem- 
perature of the star T p h is a function of 9 since energy 
balance in steady state requires 

aT^9) = aTt r {S) + F in (7) 

at each point on the stellar surface, where Fi n is the 
intrinsic energy flux coming from the stellar interior (de- 
rived from cooling of the stellar interior, gravitational 
contraction, and D burning) and 



F 



= T (Ag)^ 



(8) 



3.5 x 10 K I MtM^R^i 



1/4 
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Equation {3 assumes that all radiation intercepted by 
the star gets fully absorbed by its surface and reflection 
is negligible. Our discussion can be easily extended for 
the case of non-zero stellar albedo. 

In the absence of irradiation Fi n = uTq . With irradi- 
ation the local flux emitted by the photosphere oT^ h ex- 
ceeds cxTq , but the intrinsic stellar flux Fi n derived from 
the gravitational contraction and cooling of the stellar 
interior actually becomes smaller than ctTq as we demon- 
strate in fc|3 3.11 

We define the regime of weak irradiation as that corre- 
sponding to low M such that 

T irr (9)<T (9) 

for any 9 (i.e. Kg -C 1). As the irradiation is more 
intense near the stellar equator, weak irradiation at any 
point on the stellar surface requires Ti rr (7r/2) < Tq (or 
A < 1), or accretion rates lower than 
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Energy dissipation in the equatorial boundary layer 
(which we do not account for here) can heat equatorial 
region above To even at M < M c , but this heating does 
not spread very far from the equator and does little to 
affect the large scale stellar structure. 

A regime of strong irradiation is defined as that cor- 
responding to M > M c (A> 1) so that at least some 
parts of the stellar surface have 

T lr ,(0)>T o (11) 
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(or Kg > 1). Initially this condition is satisfied only 
near the stellar equator where an irradiated belt with 
Ti rr (6) > To forms. As M increases this belt expands in 
^-direction, although rather slowly since Fi rr is a rapidly 
decreasing function of 9, see Figure^ As will be shown in 
£ )3 3.11 in irradiated regions the intrinsic energy flux Fj„ 
coming from the stellar interior is suppressed compared 
to <tTq so that the effective temperature of the irradiated 
part of the star can be well approximated by 

T p h ~ T irr (9). (12) 
Transition between the low-latitude irradiated region and 
the high-latitude part of the stellar surface where T p h ~ 
Tq takes place at 9i„ given by (see eq. [5]) 



sm f 



RlaT£ 



1/5 



= (IiA) 



-1/5 



JxGM+M, 

« 0.5 if^R^ 5 M^ 1/5 MZl /n . (13) 

According to this formula, at M = 10~ 5 M Q yr _1 polar 
regions having temperature To ~ 3000 K occupy about 
15% of the stellar surface. The rest of the surface has 
Tph significantly modified by intense radiation coming 

from the disk. At this M equatorial temperature reaches 
T pH (-k/2) w 1.8 x 10 4 K, much higher than T - 3000 
K corresponding to the typical Hayashi track of a young 
star. 

3. COOLING OF IRRADIATED STELLAR SURFACE. 

Young stars, brown dwarfs and giant planets are fully 
convective objects. It is well known (Kippenhahn & 
Weigert 1994) that many characteristics of fully convec- 
tive objects such as their luminosity and effective temper- 
ature are determined mainly by the properties (opacity 
behavior, ratio of specific heats of the gas) of their out- 
ermost, near-photospheric layers and are rather insensi- 
tive to the processes occurring in the convective interior. 
Given that irradiation changes the boundary conditions 
on the surface of accreting fully convective object we may 
also expect that it should affect the luminosity of such 
an object (Arras & Bildsten 2006). 

Intense heating some parts of the stellar surface sup- 
presses convection in the subsurface layers and gives rise 
to a convectively stable radiative zone sandwiched be- 
tween the photosphere and convective interior, see Figure 
[1] for illustration. Appearance of this zone is analogous 
to the formation of a roughly isothermal radiative layer 
in the outer parts of the close-in extrasolar giant plan- 
ets caused by the intense radiation of their parent stars 
(Guillot et al. 1996; Burrows 2000). It is the structure of 
this zone that we want to investigate in order to assess an 
impact of irradiation on stellar cooling. Here we assume 
that the radiative zone is 

1. optically thick, as measured from its bottom 
(convective-radiative boundary) to the photo- 
sphere, and 

2. geometrically thin compared to R+. 

Validity of these assumptions is verified in £13 3.11 

In the absence of internal energy sources radiation 
transport in the optically thick radiative layer is gov- 
erned by 

V-F = 0, F = -^— VT, (14) 
3 up 



where F is the radiative flux density, n is opacity and p 
is the gas density. Equation of hydrostatic equilibrium 
reads VP = — pg, where P is the gas pressure and g 
is the local gravitational acceleration. These two equa- 
tions describe the radiative zone structure subject to the 
boundary condition 



T 



r=2/3 



T p h{0), 



(15) 



where r is the optical depth. The radiative boundary 
condition (fl5|) coupled with (|7|) is appropriate here be- 
cause most of the stellar surface is not obscured by the 
accreting gas and is free to radiate energy into space. 

3.1. ID approximation for the structure of the 
radiative layer. 

In general one must solve equation (|14p together with 
the equation of hydrostatic equilibrium in two dimen- 
sions - r and 9. However, under the circumstances clar- 
ified in £13 3.31 the r-component of the radiative flux F r is 
much larger than its ^-component Fg, so that the latitu- 
dinal transport of energy can be neglected. This leaves 
r as the only independent variable in equation (|14p ef- 
fectively making it one-dimensional. A dependence on 9 
then appears only through the external boundary condi- 
tion, namely T p h{9). This is the limit that we will focus 
on in this work. 

Because of the thinness of the radiative zone, r varies 
only weakly through the radiative zone so that one can 
neglect the divergence of the radial component of F r , 
thereby reducing equation (fT4|) to simply dF r /dr = 0. 
With these simplifications equation ([14]) can be inte- 
grated once to find 



F — — 



16 crT 3 dT 
3 np dr ' 



(16) 



where the integration constant on the left-hand side is 
independent of r and as such has to coincide with the in- 
trinsic flux Fi n coming from the convective interior of the 
star. Determination of Fi n is the goal of our calculation. 

Radial pressure gradients in the radiative zone are 
much larger than the latitudinal gradients so that the 
equation of hydrostatic equilibrium can be written as 



dP 
dr 



= ~P9, 



(17) 



where g = |g| « GM+/R+ is the gravitational accelera- 
tion which is roughly constant within the thin radiative 
layer (stellar rotation is neglected throughout this work) . 

Subsequent consideration is very similar to the calcula- 
tion of the atmospheric structure for the protoplanetary 
core immersed in a protostellar nebula, which can be 
found in Rafikov (2006). We assume that k depends on 
gas pressure and temperature as 



kP a T^ 



(18) 



where a > and fi are constants. Equation (|18|) is a rea- 
sonable approximation to the opacity behavior in some 
density and temperature intervals typical for young stars. 
In particular, at 2500 K < T < 5000 K opacity is mainly 
due to H~ absorption with electrons supplied by elements 
heavier than H with low ionization potentials. Bell & Lin 
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(1994) have demonstrated that in this regime k can be 
well fit by 



6 X l - 14 P 2 /3 T 7/ 3i 



(19) 



At T > 5000 K electrons from partial hydrogen ion- 
ization enhance H - opacity and hydrogenic absorption 
dominates. In this regime (Bell & Lin 1994) 



2.4 x i - 39 p 1 /3 T 29/3_ 



(20) 



Here we also assume that the equation of state (EOS) 
of gas in the whole star, including the external radiative 
zone, can be characterized by a single ratio of specific 
heats 7. In other words, we assume that under adiabatic 
conditions gas behaves as P = Kp 1 , where K is a con- 
stant set by the entropy of the gas and 7 is fixed through- 
out the star. We adopt 7 = 5/3, which should work fine 
in fully ionized, convective interiors of young low-mass 
stars, although this approximation is not very accurate 
at the transition between the outer radiative zone and 
the convective interior since gas is only partly atomic 
there. Continuing dissociation and ionization cause vari- 
ations of 7 in this region which may be quite important, 
see £)5 5.21 Nevertheless, to get a qualitative picture of 
the effect of irradiation on stellar cooling and for making 
rough numerical estimates this constant-7 approximation 
should be sufficient. 

With k in the form (fT8"|) equation (fTB")) can be integrated 
using (fTT)) . (fT5|) and the ideal gas law: 
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Pph 



1+a 



1 = 



ph 



T 

T p h 



4-/3 



- 1 



where 



V 



1 + a 



Vp/j — 



4-/3' 

3 Fi n Kiph Pph 

16 9°T} h 



(21) 

(22) 
(23) 



and P p h and K p h — K.Pp h T^ h are the values of pressure 
and opacity at the photosphere. 

Solution (|2"Tj) allows us to calculate temperature gra- 
dient 



1 



Tph 



4-/3 



1 



'ph 



(24) 



which determines whether gas is stable against convec- 
tion. Note that at the photosphere V(Tph) = V p /j. Ev- 
erywhere inside the radiative zone 



V < V od = (7-IV7, 



(25) 



where V at j is the adiabatic temperature gradient. In con- 
vective regions V > \7 ad . For our adopted 7 = 5/3 one 
finds V ad = 2/5. 

We also assume that at some depth an object un- 
der consideration does become convective and determine 
what is necessary for this transition to occur. If (3 < 4 
then equations (|24|) and |25|) demonstrate that convec- 
tion sets in only provided that 

V > V ad . (26) 

Situation described by equation (|2"6")) is realized e.g. for 
opacity given by equation (| 19[) . when Vo = 1 exceeds 



V a d = 2/5, implying that radiative energy transport 
does indeed change to convective at some depth, as we 
have assumed. 

On the other hand, when opacity is characterized by 
(3 > 4 equation ([24]) guarantees that transition to con- 
vection occurs at some depth, irrespective of the exact 
value of either or Vo- This situation is appropriate for 
k given by equation ([20]) since in that case (3 rs 10. 

Despite this difference, in both cases the temperature 
T c b and pressure P c b at the convective-radiative bound- 
ary are given by 4 



Pcb — Pph 



Vo - V 



ph 



Pcb = Pph [ ^ 



Vo - V ad 
Vo 



1/(4-/3) 



'ph 



Vp/j Vo — V a d 



l/(l+a) 



(27) 



(28) 



which can be derived by setting V(T c f,) = V a( 2 and using 
equation (f2"T) . Using equations (fTT)) . (|2"Tj) , and (f2"8"|) one 
can also find that the optical depth at the convective- 
radiative boundary 

-1 



T cb ' ~ v p h ■ 



V:, 1 . (29) 
while the radial extent of the outer radiative zone is 



AR r ~ H p hln\7 



ph ' 



(30) 



where H p h — ksTph/ (pg) is the photospheric scale 
height. Given that irradiation cannot heat the star to a 
temperature comparable to its central temperature (oth- 
erwise outer layers would be unbound) H p h should be 
much smaller than i?* even under rather extreme irradi- 
ation. 

Both Vo and V a d are constants of order unity. This 
makes it clear from equation (|27p that the temperature 
variation between the photosphere and the convective 



zone boundary is rather small, \T c b — T. 



ph] 



T, 



ph • 



In 



practice, we find that at T < 5000 K, when k is given by 
(fT5|) . convection sets in at T c b ~ 1.36XU, while at higher 
temperatures, when k is given by (l20j) . P c b ~ 1.19Xp/j. 
At the same time, under strong irradiation the external 
radiative zone should be deep enough for the pressure 
P c b at its bottom to greatly exceed P p h. In this case 
equation (|28f suggests that 



Vpft < 1, 



(31) 

a result that is verified in §3 3.2| see equation (|44]) . 

According to equation ([29]) smallness of V p ^ results in 
T c b S> 1, thus justifying our assumption (1) about the 
radiative zone properties. At the same time, because of 
rather weak (logarithmic) dependence of AR r on V p /i, 
the thickness of the outer radiative zone should not be 
much different from H p h -C -R*. As a result, Ai? r <C i?*, 
verifying our assumption (2). Thus, the condition (I31| 
can then be viewed as a prerequisite for the formation of 
a geometrically thin, optically thick radiative zone with 
roughly isothermal temperature profile under the action 
of intense external irradiation. External radiative layers 
with similar near-isothermal structure are expected to 
exist in the envelopes of irradiated hot Jupiters (Guillot 
et al. 1996; Baraffe et al. 2003; Chabrier et al. 2004) and 

4 These results coincide with equations (47) and (48) of Rafikov 
(2006) if we identify V p h = Voo and assume 1. 
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in the outer parts of the low-luminosity atmospheres of 
protoplanetary cores immersed into the protoplanetary 
nebulae (Rafikov 2006). 

The value of P p h can be fixed in the following way. 
Above the photosphere gas is roughly isothermal with 
temperature T p h - an approximation which is good 
enough for our purposes. Then at height z above the 
photosphere p(z) = p p h exp(— z/H p h), where p p h is the 
photospheric gas density. Using this result and equation 
(fig)) we find the photospheric optical depth 



- - f Knriz - K P' lP P h 
3 -Jnpdz- {a + i)g , 

o 



from which it follows that 



P 



ph 



2(q + 1) g 
3 < 



l/(l+a) 



(32) 



(33) 



As a byproduct of relation ()32jl one can rewrite equation 
(l23l) as 



a - 



1 Ft, 
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(34) 



that F in < 



It then follows from equations (|3"Tj) and 

We are now in position to evaluate Fi„ and see how ir- 
radiation affects cooling of convective objects. To do this 
we note that the inner boundary of the radiative zone is 
also the outer boundary of the convective interior. We 
assume that convective transport is so efficient that en- 
tropy is constant 5 throughout the inner convective zone, 
so that the EOS can be well represented by P = Kp 7 , 
where K is the adiabatic constant. As a result, P p h and 

T p h must be related via {kT c b/ p) 1 = KP^ 1 which, cou- 
pled with equations (23|), fllfj), (EH), and condition IpTTjl . 
yields the following expression for Ti„: 

V /V a 

F in (0) 



16V arf /V - V ad 



pK 1 ^ 
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(l+a)/V a 
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(35) 



where 



f = + (1 + a)/V ad . 



(36) 

Intrinsic stellar flux Fi n exhibits an explicit latitudinal 
dependence because it is a function of T p h{9). 

As discussed before, when (3 < 4 and Vo > a transi- 
tion to convection at some depth requires Vo > V Q rf. As 
a result, 

4-£=(4-/?)(l-^ <0. (37) 

On the other hand, when f3 > 4 one also finds 4 — £ < 
because Vo < in this case. Thus, in both situations Fi n 
decreases as T p h increases. In other words, irrespective 
of the opacity behavior external irradiation of the stellar 
surface suppresses stellar cooling, a result known from 

5 In reality stellar envelope contains superadiabatic regions 
which are not captured in our analysis and may affect its results. 



the studies of irradiated giant planets (Guillot et al. 1996; 
Burrows 2000). 

Since external radiative zone is rather thin compared to 
Ri, it must contain negligible amount of mass compared 
with M*. Then the structure of fully convective inner 
region of the star should be well described by the classical 
theory of polytropic spheres (Landau & Lifshitz 1984; 
Kippenhahn 1994). In particular, adiabatic constant K 
can be related to the stellar mass and radius as 



K = ^)GM^RV- 



(38) 



where C(l) ~ 1 is a parameter set by the equation of 
state of the gas. In a particular case of convective young 
stars with fully ionized interior characterized by 7 = 5/3 
one has C(5/3) = 0.1286 and 



K = 1.081 x 10 14 Ml /3 Rn. (39) 

Equations (f35|) and (|38|) unambiguously determine 
cooling of the star as a function of stellar parameters 
R± and M+, temperature distribution at the photosphere 
T p h(0), and opacity behavior in the outer radiative zone. 

3.2. Comparison with the case of an isolated star. 

We now compare stellar cooling in the irradiated case 
with that occurring in isolated stars, in the absence of 
external illumination. In the latter case T p h = To, F% n — 
F = ctTq and equation (|3"4")1 gives \7 ph = V e // = (a + 
l)/8 ~ 1. Substituting this result into equations (|27|) . 
|28p. using adiabatic relation at the convective-radiative 
boundary, and equation (|33|) we find 



T 



16V ad /V - V 



eff 



Vo - Vac 



l-Vo/V a 



9 f^K 



1/7 



(l+a)/V ad 



(40) 



This expression sets the effective temperature of the star 
and its cooling rate Fq as functions of M*, R+ and opacity 
behavior. In particular, for k typical at temperatures 
below 5000 K one finds using equation (|3T)|) that 



T 



1200 K M" /39 i?J{ 13 . 



(41) 



This is considerably smaller than T ef f sa 3000 - 4000 
K typical for an isolated fully convective star on the 
Hayashi track that one obtains with detailed numerical 
stellar structure calculations (Siess et al. 2000). We as- 
cribe this difference to our adoption of fixed 7 throughout 
the whole star and the neglect of superadiabaticity in the 
outer parts of the convective region, see £)5 5.21 At the 
same time equation (|41l) captures the main property of 
the Hayashi track - extremely weak sensitivity of To to 
-R* and, consequently, stellar luminosity. 

If we now go back to equation (|35p one can easily see 
that it can be rewritten as 



To 



T ph {9) 



or, with equation (J7]), as 
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This result together with (|37| once again vividly illus- 
trates the inhibition of stellar cooling by external irradi- 
ation and specifies the magnitude of this effect. 
Using equations ([34]) and (|42f we can also write 



V 



ph 



(a + 1) fTj 



tT 



(44) 



which shows that \7 p h <C 1 when stellar surface is 
strongly irradiated (T p h > To), thus confirming equa- 
tion (f3Tj) . Note the strong dependence of V p h on T p h/To: 
with our power-law anzatz for opacity £ « 13/2 and ~ 13 
below and above 5000 K correspondingly, see equations 
(HD and flU}. 

3.3. Conditions of validity of ID approximation. 

In Appendix [B] we determine the circumstances under 
which the results of £13 3. II hold true. We show there that 
the condition of the validity of ID approximation can be 
expressed as 



H. 



ph 



< v 



ph 3 



(45) 



where Lg is a characteristic scale in direction over which 
the external boundary condition [in our case Tgh(0)] ex- 
periences variation. Equation JT]) and Figure |2] demon- 
strate that in the case of irradiation by accretion disk 
Lg ~ R+. Then equations ([34]) and (|42]) allow us to 
rewrite the condition (|45[) as (assuming that To and T p h 
are in the same opacity regime) 



7; - 7 1 k 



2/(2+?) 



(46) 



where -ffo = fceTo/(/ig) is the photospheric scale height 
in the absence of irradiation. Given that 



^ « 5 x 10 3 AhRuT^ 



(47) 



we may conclude that ID approximation should be rather 
accurate even if T p h exceeds To by a factor of several (e.g. 
Tph % 7T for T < 5000 K). 

Whenever the condition (j4"5|) is violated the redistribu- 
tion of energy in ^-direction within the radiative layer be- 
comes important. In this case one needs to solve the full 
two-dimensional equation (| 14[) without assuming that ra- 
diative flux in 6 direction is small. A similar situation 
arises at stellar equator where a lot of energy is released 
in a boundary layer that is not very extended in 9 di- 
rection (Popham et al. 1993). As a result, at equa- 
tor Lg <C R± and the condition (|4"5)) can be violated 
there even though at all other latitudes ID approxima- 
tion works fine. 

4. INTEGRATED STELLAR COOLING. 

We are now in position to calculate the integrated in- 
trinsic luminosity L (due to stellar contraction and inte- 
rior cooling) of a convective star that is irradiated by a 
circumstellar disk: 
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L = 4nRl [ F in (9) sin 9d9, 



(48) 
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I2 computed using g(9) (see Figure [2j and the dashed line is I2 
calculated using asymptotic representation (|5j. 



where Fi n {9) is given by the expression (|35[) in the ir- 
radiated part of the stellar surface, for 9 > 6i rr , while 
Fin{9) ~ ctTq in the weakly irradiated polar regions, for 

^ Oirr- 

Convective objects can exhibit different modes of cool- 
ing which is best illustrated by considering the limit 
9irr -C 1 (A 3> 1). In this limit the contribution of 
polar caps to the total luminosity is 

(49) 
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while irradiated equatorial regions contribute 

tt/2 
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see equation (|4*2^) . 

Using equations ([5|), (fl3|) it is easy to see that for 
9irr 1 the latter integral is dominated by 9 sa 9i rr 
if 
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In this case, according to equation 
imate g{9) w / sin 5 9 and find that 



(51) 



one may approx- 



L, 
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(52) 



This results leads us to an interesting conclusion that 
as long as the condition (f5Tj) is fulfilled, an object cools 
predominantly through its polar caps and its integrated 
luminosity L is almost independent of the details of opac- 
ity behavior in its outer layers. The latter point is easy 
to understand, since in this case L ~ FqS, where S is the 
surface area of the polar caps. But according to equation 
(|13p the value of S is determined only by irradiation and 



To- As a result, L depends on k only weakly, through 
the L er contribution. 

We call the regime of stellar cooling realized under the 
condition (|5ip the high-latitude cooling. This regime nat- 
urally occurs in irradiated young stars since £ > 28/5 for 
k given by either |19|) or {20]). Equations and (|52jl 
demonstrate that in this regime L is suppressed roughly 
by ~ Q\ rr which may be as low as ~ 0.2 — 0.4 accord- 
ing to the expression (TT3]). Thus, disk irradiation can 
substantially slow down cooling of young stars. 

In the opposite case of £ < 28/5 cooling is in the low- 
latitude regime, so that star loses most of its internal en- 
ergy through the equatorial regions even though they are 
strongly irradiated. In this case one should use the full 
expression (|50[) to evaluate L w L er . Stellar luminosity 
suppression for 6 irr <C 1 is given by 

L/L ^CI 2 (sm8 zrr ) 5 ^/\ (53) 

rr/2 

I 2 (0 = lt 4)/A J[gm^ 0/i ^9d9. 
o 

Function /2(f) ~ 1 is shown in Figure S Knowing that 
4 < £ < 28/5 in the low-latitude case one can easily see 
that the degree of luminosity suppression is smaller than 
in the high-latitude cooling regime. 

In Figure |4] we plot L/Lq - the ratio of stellar lumi- 
nosities in the irradiated and isolated cases - as a func- 
tion of irradiation parameter A, for different values of £. 
This calculation does not explicitly make an assumption 
Oirr "C 1 (or A > 1), although it covers this regime as 
well. Here Lj ' Lq is computed by the straightforward in- 
tegration of Fi n over the stellar surface (including the 
polar caps where Fi n ~ Fq), with the distribution of 
Ti rr (6) found from equation ^ and g(6) displayed in 
Figure [2j We also indicate the asymptotic behavior of 
L/Lo as given by equations (j49|) and (|52]) for £ > 28/5 
and equation (|53j) for £ < 28/5. 

One can easily see that, as expected, L/L oc 6f rr oc 
A~ 2 / 5 as A ^> 1 for £ — 13 and 6.5 independent of the ac- 
tual value of £ (only the normalizations of the curves are 
different because of the different contributions produced 
by the near-polar cap regions) since for both £ > 28/5. 
Significant suppression of the stellar flux (by ~ 2) is 
found in this case already at A ~ 10 2 — 10 3 . In the case 
£ < 28/5 asymptotic behavior for A> 1 agrees well with 
equation (|53|) . L/L oc A^ 4- ^/ 4 , and the degree of stel- 
lar flux suppression is weaker than in the high-latitude 
regime: L/L ~ 0.5 only at A » 5 x 10 3 for £ = 5.3 and 
at A « 2 x 10 5 for £ = 4.5. 

Note that results presented in Figure [3] are calculated 
neglecting any additional heating that can be produced 
near the equator by the boundary layer dissipation. We 
will address this point in more detail in ^ 35 5.21 

5. DISCUSSION. 

Luminosity suppression by disk irradiation is one of 
the most important results of this work. Analogous phe- 
nomenon has been previously found in studies of extra- 
solar giant planets in short-period orbits, where stellar 
irradiation is quite severe (Guillot et al. 1996; Burrows 
et al. 2000). In that case irradiation affects only one 
side of the plant which always faces the star. Some heat 



iiiiiii| iiiiiiii| lining iiiiiiii| iiiiiiii| iiiiiii^ iiiiiiii| iiiiiiij iiiiiii^ iiiiiiii| iiiiiii| iiiiiiij 111111 




1 10 10 s 10 3 10 4 10 6 10 6 10 7 10" lOMO'oiO 11 



A=GM.M/(crT 4 , R?) 

Fig. 4. — Plots of the intrinsic luminosity of an irradiated star 
L (in units of an isolated star luminosity Lq) as a function of irra- 
diation parameter A = GM+M / (aT^ R +) fo r different values of the 
power law index £ defined by equation l|36|l . Dotted lines illustrate 
the corresponding asymptotic behaviors for A ^> 1: L/Lq oc A -2 / 5 
for £ > 28/5 and L/L oc A( 4 ~^/ 4 for £ < 28/5. 



from the day side gets redistributed to the night side 
by atmospheric circulation (Menou et al. 2003; Dobbs- 
Dixon & Lin 2007) which should complicate the calcula- 
tion of the photospheric boundary conditions across the 
whole planetary surface. In our case irradiation is az- 
imuthally symmetric which makes our calculation more 
robust. Analogous to the case of extrasolar giant plan- 
ets we expect that luminosity suppression by irradiation 
would tend to retain heat inside the star and increase 
stellar radius above the value found in the absence of 
irradiation (Baraffe et al. 2003). Whether this radius 
increase is significant will be investigated in future work. 

Luminosity suppression may have some effect on the 
strength of the magnetic field that is generated by dy- 
namo action in the convective interior of the star. Since 
in irradiated case convective eddies transport smaller en- 
ergy flux to the stellar surface than in the case of an 
isolated star the speed of convective motions is expected 
to be smaller. This results in a less vigorous dynamo 
action and likely weaker magnetic field generated inside 
the star. 

Formation of an optically thick radiative zone near the 
stellar surface in irradiated regions is another result of 
this work which has important implications. As the star 
grows the hot gas in the vicinity of the boundary layer 
where disk meets the star gets advected into the con- 
vective interior thereby raising stellar entropy. This pro- 
vides another way of slowing down stellar contraction, in 
addition to the luminosity suppression discussed above. 
Hartmann et al. (1997) argued that impact of the heat 
advection on stellar structure is not significant as long as 
the temperature of advected gas is much smaller than the 
central temperature of the star. In the irradiated case 
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it is the properties of the external radiative zone that 
determine the temperature of the gas at the convective- 
radiative boundary, and thus the amount of thermal en- 
ergy advected into the stellar interior. Indeed, hot gas 
sinking through the radiative zone will lose a significant 
fraction of its thermal energy by radiative diffusion in 
the latitudinal direction, so that the temperature at the 
convective-radiative boundary is likely to be lower than 
in the center of the boundary layer. It is important to 
build the detailed 2D model of the radiative transport in 
the vicinity of the boundary layer to quantify this effect 
and to verify the significance of heat advection (see also 

Presence of the radiative zone may also affect atmo- 
spheric opacity in accreting brown dwarfs and giant plan- 
ets. Accreted gas brings in significant amount of dust 
into the object's atmosphere (dust can also form out of 
the gas phase under low-temperature conditions) which 
changes k and radiative properties of the star (Chabrier 
et al. 2000). However, if the object is fully convective all 
the way to its photosphere, vertical fluid motions quickly 
advect dust grains into the hot interior where grains get 
easily destroyed. This is not the case in irradiated re- 
gions of accreting objects since dust grains can hover in 
the radiative zone for a long time, as long as their grav- 
itational settling is not too fast. Of course, for grains to 
exist in the outer radiative layer in the first place Ti rr 
must be lower than the sublimation temperature of dust 
material, which may be possible only in accreting brown 
dwarfs and giant planets. 

We expect that in the case of young stars it would be 
very difficult to obtain a direct observational confirma- 
tion of the luminosity suppression by irradiation. The 
major reason for this is that L/Lq starts to deviate from 
unity only when M and, correspondingly, accretion lu- 
minosity are very large. At this stage the luminosity of a 
star+disk system is completely dominated by the direct 
emission from the disk and the disk flux intercepted and 
re-radiated by the stellar surface. Intrinsic stellar lumi- 
nosity provides negligible contribution which would be 
almost impossible to distinguish. Besides, forming pro- 
tostar should still be enshrouded in the dense veil of the 
residual gas collapsing onto the circumstellar disk. Re- 
processing of star+disk emission in this infalling envelope 
would complicate things even more. Another potential 
way of detecting the luminosity suppression is indirect, 
through its effect on the stellar radius and luminosity as 
the star emerges as almost fully formed Class I object in 
the end of active accretion phase. 

Effects of disk accretion in star formation have been 
previously investigated by a number of authors. Adams 
& Shu (1986) and Popham (1997) have calculated the 
amount of energy which is emitted by the disk and is 
intercepted by the star. Unlike us these authors were not 
primarily concerned in the details of the distribution of 
irradiation flux over the stellar surface. This is a crucial 
point of our study allowing us to identify the two different 
regimes of stellar cooling - high- and low-latitude. 

Mercer-Smith et al. (1984) were the first to explore 
the effect of disk accretion on the stellar structure. They 
handled disk accretion by specifying mass addition rate 
and accretion luminosity as external boundary condi- 
tions. They find that the stars formed by disk accre- 



tion have larger radii than nonaccreting stars of the same 
mass. As mentioned in Hartmann et al. (1997) this out- 
come most likely results from allowing the accreted ma- 
terial to have very high entropy which leads to stellar 
swelling, see Prialnik & Livio (1985). This approxima- 
tion is unlikely to be valid in reality since disk material 
joining the stellar surface should have enough time to ra- 
diate most of its thermal energy before being fully incor- 
porated into the star. Palla & Stahler (1992) and Hart- 
mann et al. (1997) in their studies of intermediate- and 
low-mass stars allowed the accreted material to have low 
entropy. They found that accretion reduces stellar size 
compared to the non-accreting case since in this case the 
addition of mass leads only to the increase of gravita- 
tional energy of the star and is not accompanied by the 
increase of thermal energy. 

All these studies have either ignored irradiation of star 
by the disk or accounted for it only in the averaged sense 
which may not be acceptable as our study demonstrates. 
To get a complete picture of protostellar evolution one 
needs to include the luminosity suppression by disk irra- 
diation into account. Such calculation must necessarily 
allow for the spatial distribution of irradiation flux on the 
stellar surface since only in this way a proper estimate of 
the luminosity suppression can be obtained. 

5.1. Applications to real systems. 

Here we apply our results to different classes of fully 
convective objects which may accrete through the disk at 
high M . In doing our estimates, which require the knowl- 
edge of R+ and T , we use the radii and photospheric 
temperatures of corresponding objects determined in the 
absence of irradiation and mass inflow by accretion as 
proxies for i?* and To that these objects would have if 
irradiation and accretion were properly accounted for. 
Needless to say, a truly accurate estimate of the effect of 
irradiation can be obtained only if i?* and Tq are calcu- 
lated self-consistently accounting for the effects of accre- 
tion and irradiation. 

5.1.1. Young stars 

Young low-mass stars transitioning from Class to 
Class I phase are fully convective and should be assem- 
bled by mass accretion from a circumstellar disk within 
several 10 5 yrs. This implies very high accretion rate 
and we adopt M = 5 x 10~ 6 M© yr _1 for a simple es- 
timate. An isolated M+ = 0.5 M star at an age of 10 5 
yrs has a radius R+ — 3.9 R and effective temperature 
T = 3760 K (Siess et al. 2000). These parameters yield 
A = 10 2 and T lrr (n/2) w 6000 K, although the latter 
is likely to be higher because of the boundary layer dis- 
sipation. At this A cooling in the equatorial region is 
suppressed but the size of the polar caps is reduced only 
weakly: according to Figure is given 

by an implicit relation g(6i rr ) = A -1 ]. In Figure we 
present more general results for A calculated for stars of 
different masses using stellar parameters from Siess et al. 
(2000) and assuming constant M = M*/i acc , where t acc 
is the accretion time (assumed equal to the stellar age). 
One can see that the low-mass stars assembled within 
3 x 10 5 yrs generally have A in the range of 30 — 10 2 , 
agreeing with our simple estimate. 

Since £ w 6 for T < 5000 K, young stars cool mainly 
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through the polar regions and we find from Figure [4] that 
stellar luminosity is reduced by irradiation only by about 
10% for A sa 10 2 . On the other hand, if M is not con- 
stant but increases as M* grows one may expect values 
of A larger by a factor of several. Also, to calculate A we 
have used parameters of isolated stars while stars assem- 
bled by disk accretion of the low-entropy material have 
smaller R+, leading to larger A. All these factors may in- 
crease the importance of disk irradiation in determining 
L of young accreting stars. 

5.1.2. Young stars in quasar disks 

A very interesting mode of star formation is possible in 
the accretion disks around the supermassive black holes 
(SMBHs) in the centers of galaxies (Illarionov & Ro- 
manova 1988; Goodman & Tan 2004; Nayakshin 2006). 
It is currently known that our own Galactic Center har- 
boring a SMBH of mass M BH w 3.7 x 10 6 M© (Ghez ct 
al. 2005) contains a number of young (< 6 Myrs) mas- 
sive (M* > 10 Mq) stars that form two misaligned disks 
around the SMBH (Paumard et al. 2006). One of the 
most likely scenarios for the origin of these stars is a 
fragmentation of a gravitationally unstable gaseous disk 
(or disks) followed by the growth of fragments to their 
present masses by gas accretion from the residual disk 
(Levin 2006; Nayakshin 2006). Assuming that disk tem- 
perature is kept at the level of 50 K by the radiation of 
nearby stars (Levin 2007) one finds that at a = 0.1 pc 
from the SMBH (which is the typical dimension of the 
observed disks) surface mass density of £ sa 27 g cm~ 2 
is required for the disk to be Toomre unstable. 

Fragments formed as a result of instability at 0.1 
pc have a typical mass M+ ~ S/i 2 w 10 -3 M Q (ap- 
proximately one Jupiter mass) where ft, is a disk scale 
height. At formation the Hill radius of such an ob- 
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Fig. 5. — Irradiation parameter A for stars of different mass and 
age t aC c (indicated on the plot) assembled by disk accretion with 
constant M = M+ftacc- Stellar parameters from Siess et al. (2000) 
were used in this calculation. 



ject Rh = a(M+/MBH) 1 ^ 3 is already comparable to 
h and as M* grows by accretion Rh becomes larger 
than h. As a result, accretion onto fragment proceeds 
through the sub-disk that forms within the fragment's 
Hill sphere, presenting us with the setting investigated 
in this paper. Rate at which gas flows into the frag- 
ment's Hill sphere is the Hill accretion rate 6 Mg « 
mR% ~ 2 x 10~ 4 Mf /3 (a/0.1 pc) 2 M yr" 1 . Note 
that Mh is smaller than the Eddington mass accretion 
rate M E dd = kncR+l 'n es = 1.4 x 10~ 3 i?n M yr -1 (here 
c is the speed of light and K es is the electron scattering 
opacity) at a = 0.1 pc but may become comparable to 
MEdd further out from the SMBH provided that R+ is 
not much larger than R@ . 

If gas in the disk is able to accrete at the same high 
rate Mh onto the stellar surface then 

A w 3 x 10 5 M 1 5/3 T 3 7 5 4 i? n 3 . (54) 

At present we do not have a theory for the structure 
of stars formed by fragmentation of gravitationally un- 
stable disks so the value of i?* is highly uncertain. If 
R* < 10 Rq then A > 10 3 and luminosity suppression 
by irradiation should be quite important, reducing L by a 
factor of 2 — 3 compared to Lq, as Figure [4] demonstrates 
for £ = 6.5. As the value of i?+ itself is affected by the 
time history of L, irradiation should not be overlooked 
in studies of star formation in quasar disks. 

5.1.3. Young brown dwarfs 

Brown dwarf (BD) formation is likely to be a scaled 
down version of the low-mass star formation: one again 
expects a formation of a centrifugally supported disk 
around a fully convective object that grows by disk accre- 
tion. The biggest uncertainty in determining A is again 
i?*: 0.1 Gyr old BDs have radii of 0.1 - 0.2 R (Baraffe 
et al. 2003) but accumulation of their mass (poorly in- 
vestigated at present) likely takes less than 10 5 yr, dur- 
ing which time their entropy is still quite high, result- 
ing in considerably larger R+. Assuming that an ob- 
ject with M+ = 0.03 M Q grows at constant M in time 
t acc = 5 x 10 4 yr and has R+ = 0.5 R Q and T » 3000 K 
we find A w 800. Provided that opacity can still be char- 
acterized by expression (|19[) (which is a somewhat ques- 
tionable assumption) we conclude that irradiation may 
lead to order unity reduction in L. As the brown dwarf 
cools and contracts A increases bringing down L/L even 
more, provided that M could still be maintained at high 
level. Thus, young BDs may be affected by the disk 
irradiation which may have consequences for their sub- 
sequent thermal evolution. 

5.1.4. Young giant planets 

Finally, we consider the situation arising during the 
late stages of giant planet formation via the so-called 
core instability. This scenario of planet formation as- 
sumes buildup of a ~ 10 M refractory core in the pro- 
toplanetary nebula by planetesimal agglomeration. The 

6 Such high M is also typical for FU Orioni objects. As demon- 
strated by Popham et al. (1993) in the high-M regime the boundary 
layer is so thick that it covers a significant (~ 0.5) fraction of the 
stellar surface slowing down interior cooling. Heat advection into 
the stellar interior may also become an issue (Popham 1997). 
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self-gravity of the core triggers an instability and leads to 
rapid gas accumulation (Mizuno 1980). While the initial 
stages of this process can be adequately described in the 
spherically-symmetric approximation, the later epoch of 
unstable gas accretion must have distinctly non-spherical 
morphology. Indeed, as mentioned in Rafikov (2006), as 
soon as the mass of a rapidly growing planet exceeds the 

so-called transitional mass M tr = c 3 /f2G w 40M^a^ 4 
(here c s is the gas sound speed in the nebula and as = 
a/5 AU is the planetary semi-major axis scaled by 5 AU) 
the Hill radius of the planet Rh becomes larger than the 
scale height of the disk h. As a result, protoplanet starts 
accreting gas from the surrounding nebula through the 
sub-disk that forms within its Hill sphere, thereby pre- 
senting a situation analogous to the star formation in 
the Galactic Center described in £15 5.1 5.1.21 (except that 
now the collapsing fragment of a gravitationally unsta- 
ble disk is replaced by a growing planet). Here we assess 
how important can irradiation be for planetary cooling 
when M p > M tr . 

The maximum M available to the planet is still likely 
given by the Hill rate M H ~ 2 x lO^AT^fa^ 1 Mj yr" 1 , 
where M p ^ = M p /10 2 M ffi and we have adopted a sur- 
face density profile S = 270a 5 3 ^ 2 g cm" 2 typical for the 
Minimum-Mass Solar Nebula. This allows us to compute 

^ A ~ 104 <V (A)" 3 (55, 

This estimate is rather uncertain because of poorly con- 
strained R p and To during the stage of active gas ac- 
cretion by the planet. Here we adopt R p — 5 Rj and 
To = 10 3 K mainly for illustrative purposes. 

Dust is the major source of opacity in the outer layers 
of forming giant planets. Dust opacity scales as k cx T 13 
with P w 0.5 — 2 depending on dust grain composition, 
spectrum of grain sizes, etc. Here we adopt j3 = 1 in 
which case £ = 4.5. This corresponds to cooling domi- 
nated by the equatorial regions, which is different from 
the stellar case. As a result, L/Lq should be more sensi- 
tive to the structure of the boundary layer through which 
disk material accretes onto the planet, namely, L/Lq 
should be lower than Figure [55] implies. Forgetting about 
this complication for the moment and using £ = 4.5 and 
A from equation (|55|) . we find from Figure [4] that in the 
planetary case T can be suppressed compared to To by 
several tens of per cent. Given the existing uncertain- 
ties in modeling the late stages of planet formation this 
degree of luminosity suppression by irradiation may not 
seem like a serious issue. However, more massive, com- 
pact and cooler planets can easily have A ~ 10 5 — 10 6 
in which case irradiation would reduce T by a factor of 
several potentially affecting planetary evolution. 

5.2. Additional complications. 

Here we address various complications that may arise 
when the results of this work are applied to real objects. 

All our derivations and estimates explicitly assumed 
opacity in the form (fT8|) . While this representation can 
be quite accurate within some temperature intervals one 
has to bear in mind that on the surface of a star ir- 
radiated by the disk temperature can vary appreciably 
between the equator and the poles. Indeed, the equato- 
rial regions of a AT = 1 M Q , R+ — 2 R Q star accreting 



at M = 10" 5 Mq yr" 1 are heated to 1.1 x 10 4 K so 
that equation ([20]) applies, while polar caps still have 
T w 3 x 10 3 K so that equation ([TO]) is more appropriate. 
Thus, in different regions of stellar surface k has differ- 
ent dependence on P and T. In this case equation (153|) 
becomes invalid and to properly compute the luminos- 
ity suppression one would need to take into account the 
latitudinal variation of not only T rr but also k(P,T). 

Even at a fixed latitude opacity can switch from one 
regime to another within the outer radiative zone. Al- 
though k is much more sensitive to T than to P and the 
radiative zone is roughly isothermal, pressure at its bot- 
tom P c b is ~ 1+ ^ 1 of the photospheric pressure 
Pph, so that even a weak dependence of k on P can lead 
to opacity transition within the outer radiative zone. In 
particular, this situation is likely to occur at latitudes 
where T irr w 5000 K and n switches from ([19]) to ([20]) . 
In this case calculation of the local radiative flux Fj„ 
gets more complicated as the external radiative region 
splits into two layers characterized by different opacity 
behaviors. 

Our numerical estimates of the luminosity suppression 
strongly rely on the assumption of fixed 7 within the 
radiative zone. Equation (|4ip which assumes 7 = 5/3 
throughout the whole star fails to predict the correct 
photospheric temperature of an isolated star. The reason 
for this is the variation of 7 at low T caused by molec- 
ular dissociation. This results in V ' a d 0.1 and leads 
to a smaller drop of temperature in the outer convective 
parts of the star. On the other hand, superadiabaticity 
of convection in the outer layers of convective zone coun- 
teracts this effect to some extent. Our simple estimate 
([4Tj) underpredicts To by a factor of 2 — 3 which is pri- 
marily a consequence of our assumption of fixed 7 = 5/3. 
A proper calculation of Fi n (0) and L/Lq must be able to 
account for the variation of 7 with P and T inside the 
outer radiative zone. 

Our analysis is affected to some extent by the presence 
of the boundary layer through which disk material joins 
the star. Viscous dissipation in this layer heats accreting 
gas to very high temperature. Since this energy release 
takes place very close to the stellar surface most of the 
heat is likely to leak out and not get carried into the 
star with the accreted gas. However, some residual heat 
may still be accreted. Moreover, in addition to advective 
there could also be a radiative energy transfer from the 
boundary layer into the outer layers of the star (Popham 
1997). Increase of T in the external radiative zone driven 
by these processes acts to additionally slow down stellar 
cooling in the equatorial region, as equation (|35p demon- 
strates. Thus, presence of the boundary layer reduces 
L/Lq even more than our analysis predicts, and the re- 
sults presented in Figure 2] should be viewed as upper 
limits on L/Lq. This effect is likely not very important 
for young stars which cool predominantly through their 
polar regions, largely unaffected by the additional heat 
deposition at the equator. However, in the case of young 
giant planets which lose fair amount of energy through 
the low-latitude part of the surface (see §5 5.1 5.1 .4|) the 
reduction of intrinsic flux in the equatorial region may 
produce quite noticeable decrease of L/Lq compared to 
the idealized case considered in this work. 

When discussing the external radiative zone through- 
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Fig. 6. — Schematic representation of magnetically channeled 
accretion. Disk material travels along the magnetic field lines and 
sediments onto the magnetospheric columns (shaded) heated by 
shock dissipation and gravitational settling of material in them. 
Hot magnetospheric columns irradiate stellar surface near the mag- 
netic poles. 

out this work we have been concerned only with the ra- 
diative energy transport. At the same time, it is well 
known that in the case of hot Jupiters advective trans- 
port in the form of atmospheric jets and winds can be 
quite important in redistributing heat across the plan- 
etary surface (Menou et al. 2003; Dobbs-Dixon & Lin 
2007). In our azimuthally-symmetric setting only merid- 
ional atmospheric motions can lead to energy exchange 
between the hot equatorial and cold polar regions. Fluid 
motions occurring on surfaces of constant effective po- 
tential (gravitational plus centrifugal) are unlikely to pro- 
duce efficient equator-pole energy exchange (compared to 
the radiative transfer) because of rather fast rotation typ- 
ical for objects formed by disk accretion. Rotation forces 
angular momentum conservation and prevents significant 
fluid motions in ^-direction, suppressing this mode of ad- 
vective transport. On the other hand, rotation tends to 
promote meridional circulation within the radiative layer 
(Kippenhahn & Weigert 1994) whose impact on the en- 
ergy transport in the outer radiative zone should be in- 
vestigated in more detail. 

Finally, our basic assumption of direct mass accretion 
from the disk may be challenged if the growing star pos- 
sesses magnetic field strong enough to disrupt accretion 
flow outside of i?* (Konigl 1991; Matt & Pudritz 2005). 
In this case gas is channeled by the magnetic field and 
is deposition onto the stellar surface at magnetic poles 
rather than at stellar equator. This completely changes 
the topology of accretion flow but the major conclusions 
about the effect on stellar cooling are likely to hold. In- 
deed, the magnetically channeled gas travels towards the 
stellar surface at a good fraction of the free-fall veloc- 
ity and at some point it must pass through the radiative 
shock, after which it accumulates at the top of the mag- 
netospheric column of accreted material, as schematically 
indicated in Figure [BJ Total energy release within the 
shock and magnetospheric column is comparable to that 



occurring if the accretion disk were extending all the way 
to the stellar surface. This hot column of accreted mate- 
rial illuminates the surface of the star leading to the same 
suppression of intrinsic stellar flux as we discussed in this 
work. In this case, however, irradiation is strongest near 
the magnetic poles while the magnetic equator is likely 
to be the coolest part of the stellar surface 7 . Calculation 
of stellar irradiation and integrated luminosity in this 
case would involve constructing a model for the magne- 
tospheric column structure and its radiative properties. 

The impact of these details on the structure and evolu- 
tion of young stars should be addressed by future work. 

6. SUMMARY. 

Luminosity of young stars actively accreting from the 
circumstellar disk can be significantly affected by the ra- 
diation which is produced in the inner parts of the disk 
and is intercepted by the stellar surface. We showed that 
if a star gains its mass via disk accretion on timescale 
of several 10 5 yr then the radiative flux caused by vis- 
cous dissipation in the disk is more than sufficient to 
increase the surface temperature of the star above the 
photospheric temperature that an isolated star with the 
same mass and radius would have. Irradiation by the 
disk is strongest in the equatorial regions and is almost 
negligible near the poles. An outer radiative zone of al- 
most constant temperature forms above the fully convec- 
tive interior in the strongly irradiated parts of the stellar 
surface. This leads to the local suppression of intrinsic 
energy flux escaping from the stellar interior. 

We have demonstrated that there are two distinct 
modes in which a fully convective object can cool: mainly 
through the cool high-latitude polar regions or predom- 
inantly through the low-latitude parts of the stellar sur- 
face. A particular regime of cooling in a given object is 
set by the opacity behavior and the adiabatic tempera- 
ture gradient V Qt ; in the outer radiative zone. Accreting 
young stars and brown dwarfs cool mainly through the 
polar regions while forming giant planets cool through 
the whole surface. 

Integrated stellar luminosity in accreting case is sup- 
pressed compared to the case of an isolated object, by 
up to a factor of several in some classes of objects (ac- 
tively accreting brown dwarfs and planets, stars forming 
in gravitationally unstable disks in galactic nuclei) . This 
leads to larger radii of irradiated objects and may af- 
fect the initial conditions which are used to calculate the 
evolution of the low-mass objects on timescales of ~ 10 
Myr after their formation. Existence of external radia- 
tive zone may facilitate retention of dust in the atmo- 
spheres of brown dwarfs and planets, and may affect the 
strength of magnetic field generated by internal dynamo 
in convective objects. 

Some of the results obtained in this work may be appli- 
cable to accreting white dwarf and neutron star systems. 

I am grateful to Gilles Chabrier for careful reading of 
the manuscript and many useful suggestion. The finan- 
cial support for this work is provided by the Canada Re- 
search Chairs program and a NSERC Discovery grant. 
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7 Illumination of the star by the distant parts of accretion disk, with the distance from the star 

beyond the point where accretion flow is disrupted by the magnetic 
field, is unlikely to be very important given the rapid fall-off of 

APPENDIX 
IRRADIATION FLUX. 

Performing an integral over </> in eq. ((T|) we find 

oo 

2 cos 9 f Fd(xR+)q(x,9)xdx 



F irr (9) = I ^ ' ' , (Al) 



/ m / 9 ■ 2/) — 7 n^C 1 - 2sin 2 6») + 1 / x 2 + 2xsin0 + 1 / x 2 sin 2 9 - i i 

qi x,9) = V&^I-2^ ^^^arctan ^ gsinfl + 1 ^ D , 1 , (A2. 

where D 2 = (x 2 + l) 2 — Ax 2 sin 2 6* and Xi n = 1/ sin.0. For Fd(R) obeying (JSJ) equation (lAlj) can be rewritten as eq. ((3j) 
with 

oo 

3 cos 6» /" f{xR i ,)q{x,9)dx 
9{6) = 1^" J x^D 2 (A3) 

Integral in (|A3[) is dominated by x ~ Near equator, where 9 — > 7r/2 one can expand integrand in (|A1[) in terms 
of 7r/2 — -C 1 and i - 1 « 1 which results in Fd(Tv/2) — Fd(R*)/2. Polar regions of the star (9 — > 0) arc illuminated 
only by distant parts of the disk, R > R+/ sm.9 ^> i?*, so that x ^> 1 (while xs'm9 ~ 1) in equation (|A1|) . Also, far 
from the star one can safely use equation |(5J| with / = 1 to finally arrive at the equation J5]) with 



h = — o / -g- \A 2 - 1 - 2 arctan = ^ (A4) 

4tt 2 7 i 6 ^ V 1 + 1 J 50tt 2 V ' 

This result is independent of the structure of the boundary layer near the stellar surface since the polar regions of the 
star do not have direct sight lines to the boundary layer. This is evidenced by the convergence at 9 — > of the two 
curves in Figure [2] calculated assuming f(R) = 1 and f(R) = 1 — (R+/R) 2 . 

VALIDITY OF lD APPROXIMATION. 

To determine the validity limits of the ID solution for the structure of the radiative zone found in £13 3.1l we evaluate 
the magnitude of the corrections arising when the latitudinal radiative transfer is accounted for. Considering ID 
solution (|2ip as a zeroth-order approximation we plug it into the full equation (|14[) and carefully expand all 9- 
derivatives, remembering that P is almost independent of 9 (latitudinal pressure gradients are small). Integrating the 
resultant expression once over r we again arrive at the equation (|16[) but with Fi n — > Fi n + 5Fi n in the left-hand side, 
where 




W*» = —7^ -5^Z(P)dP. (Bl) 



Here Z(P) ~ 1 is a weak function of pressure (varying by at most a factor ~ 1) and Lg is a characteristic scale of 
latitudinal variation of T p h, Lg — R ir (d\nT p f l /d9)~ 1 . Our ID approximation is justified if the correction to the ID 
result 5F in is small compared to F in given by equation (|35[) . 

Integral in (|B1[) attains its highest value at P ~ P p h (latitudinal radiation transfer is easiest in the upper, low density 
layers of the star just below the photosphere) and one can easily find using equations (f33|) . f|34f). and (f35|) that 

dF '" fHpH ^V-l (B2) 



Fin \ Lg J ph 

where H p h is a photospheric scale height. This result makes it clear that the ID solution for the structure of the 
radiative zone should be reasonable as long as the condition (|45|) is fulfilled. 
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